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The computation of multifractal scaling properties associated with a critical field theory involves 
non-local operators and remains an open problem using conventional techniques of field theory. We 
propose a new description of Gaussian field theories in terms of random Cantor sets and show how 
universal multifractal scaling exponents can be calculated. We use this approach to characterize the 
multifractal critical wave function of Dirac fermions interacting with a random vector potential in 
two spatial dimensions. 
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Studies of critical theories are usually concentrated 
on scaling properties of local operators and their corre- 
lations. However, there are objects in critical systems 
(at least in critical theories of disordered systems) that 
have very complex scaling properties, namely multifrac- 
tal They are characterized not by a single scaling 
exponent, but by infinitely many scaling dimensions. An 
open question of field theory is how to describe such com- 
plex scaling behavior, in particular how to calculate the 
multifractal scaling functions r(q) or f{a) 

Multifractality is a property of a distribution in a crit- 
ical theory or some other deterministic systems, such as 
dynamical systems [^1). The key feature of the phe- 
nomenon of multifractality is the property of selfsimi- 
larity of the distribution at different length scales. A 
more refined description of the property of selfsimilar- 
ity relies on the concept of box distribution. Let O be 
a large box of volume L'^, and let ijl{x) be a normal- 
ized distribution (assumed to be strictly positive) de- 
fined on 17. Imagine dividing 17 into M boxes 12^, each 
of volume a'*. The box distribution is then defined by 
= df'x lJ,{x), i — 1, . . . , M}, and the distribution 

/i(x) is selfsimilar if the scaling law 



Tiq)^D{q)iq~l) 




holds for real q ||]. The function T{q) can be shown to be 
concave and strictly increasing ||^. If D{q) is constant, 
the distribution /i(a;) is said to be simple fractal of D 
(Hausdorff) dimensional support, li D{q) is not constant 
[r(q) nonlinear], the distribution fJ.{x) is multifractal. 

The scaling function r(g) is known for many examples 
of deterministic multifractals (|,^. However, much less 
is known for distributions in critical statistical systems. 
Statistical multifractals can be found in the problem of 
localization in two spatial dimensions where the distribu- 
tion /x(a;, V) is constructed from the squared amplitude 
of the wave function solving Schrodinger equation with 
an impurity potential V{x) In this context, a very 
important example of multifractality is that of the criti- 
cal wave function at the plateau transition in the integer 



quantum Hall effect (IQHE). Numerical calculations of 
T;q[jj,(g) for the IQHE have been performed for different 
realizations of the impurity potential. Within numerical 
errors, T^qYicio) is independent of the disorder More- 
over, the exponent Tjqjj^(2), which governs anomalous dif- 
fusion 1 10 1, is consistent with the value ly fa 7/3 for the 
localization length exponent [0. Very little is known 
on the analytical structure of the critical theory in the 
IQHE. Much more is known, however, on that for Dirac 
fermions in two spatial dimensions interacting with the 
random vector potential A^{x) = e^^pOpQ^x) -\- d^x{^)i 
distributed according to exp[— 2^ J(PxA^]. In this prob- 
lem of localization, one can construct a distribution 



^{x, 8) oc 



, P[e]cxe '"A 



(2) 



from the zero-energy eigenfunctions ^{^{x) cx exp[0(x) -I- 
ix{x)\ Note that in Eq. (||) we are dealing with 

a statistical ensemble of distributions with weights ^'[O]. 
Although the disordered critical point is exactly solvable 
|]l2| , we still lack a valid calculation of Tqjj.^^((7) for the 
distributions in Eq. (||). This is so due to the failure 
of the crucial assumption commonly used in field theo- 
retical calculations of multifractal scaling exponents , 
(j^ whereby it is assumed that averaging e^'®'^' over space 



is equivalent to averaging over the disorder ||6|,|12|. Un- 
der this assumption, one may replace the inverse of the 
normalization factor J d^xe^^'^^ by a number, and use 
conventional methods from field theory. But, the result- 



ing 



{q) does not satisfy the analytical properties 



of Eq. (|1|). When e^®^^^ is not self averaging, conven- 
tional methods can not treat random variables such as 
[/ d^xe^®(^)]~^ anymore. 

One goal of this letter is to calculate T^ij^^{q) exactly 
for the first time. We hope that by doing so, we can gain 
some familiarity with the unusual properties of critical 
phenomenon for disordered system, in particular those 
associated with the plateau transition in the IQHE. Since 
conventional methods relying on the underlying critical 
field theory are unsatisfactory [p^ , we take a very differ- 
ent approach. We first construct explicitly a model for a 
statistical ensemble of distributions /x(x, ($ labels the 



1 



members of the ensemble) and then show that this model 
is closely related to the critical theory of Dirac fermions 
with random vector potentials as described, say, in . 
Our model is a simple generalization of the deterministic 
two scale Cantor set multifractal [H^. Hence, we call it 
a random Cantor set (RCS) construction. It should be 
said that the model itself has been used, for example, in 
the study of turbulences j^] and of directed polymers in a 
random medium . The novelty of our work is to relate 
the RCS to the Gaussian field theory L = i(V$)^, and 
then use this relationship to calculate T-Q-^^^iq)- In this 
way, we hope to illustrate how a new type of equivalence 
between critical field theories can be taken advantage of. 

Our first result is to show that a particular RCS con- 
struction describes a Gaussian field theory on an ultra- 
metric space. The correlation functions share the same 
form as that in the theory C — ^(V^)^ with an Eu- 
clidean metric. Although the two metrics differ, scaling 
exponents of correlation functions are the same in the two 
theories. Thus, we can use the one realization in which 
the calculation is simpler. In the case of the multifrac- 
tal scaling exponents T{q), the RCS is the realization of 
choice. 

Using the RCS realization, we can then show that 
almost all members <&) of the statistical ensemble 
yield, with the help of Eq. (|l|), a scahng exponent r(g, $) 
which is self averaging, i.e., is almost surely independent 
of $ in the limit ^ 0. To derive this result, we borrow 
an exact solution for the free energy of directed poly- 
mers on a random tree showing that the free energy is 
selfaveraging [ p^|jT^ ]. 

We provide support for the claim that the scaling ex- 
ponents r(q) should be the same in both theories by com- 
paring exact results from the RCS construction with nu- 
merical simulations for the C — ^(V$)^ theory. We also 
present analytical arguments which indicate a phase tran- 
sition in T{q) for the C — ^(V^)^ theory at exactly the 
same q^ for which a transition occurs in the RCS model. 

We begin with the construction of the RCS model. 
Consider a box Q, of volume in two dimensional Eu- 
clidean space. Imagine dividing Q. into two boxes of equal 
volume and assign to each half a binary address 

Si = ±, or, equivalently, the address (1,1) and (1,2). 
The procedure is then iterated n times. At level n the 
space is divided into 2" small boxes addressed by either n 
bits {si, . . . , s„) taking the values Sj = ±, j = 1, . . . , n or 
by (n, j) where j = 1, . . . , 2". Any arbitrary point a; G 57 
is then uniquely addressed according to the box in which 
it hes, say, by the binary address (sf , . . . , s^). Having 
constructed a selfsimilar structure, namely a binary tree, 
we define on it random variables. At each level or genera- 
tion i = 1, . . . , n of the binary tree one draws 2* indepen- 



dent random variables 



from a Gaussian 

^1/2, 



probability distribution = ( 5^)^^^ 6^P(~ fg ) (^^^ 

Figure |l|). Finally, for any given binary tree made of n 
generations, we define the random variable 



'i>„(^) = E' 



E 



2' 



(3) 



Here, for any generation i there exists one and only one 
integer < ii{x,,) < 2* such that e^jix,.) = ^jj^ix^)- 
By construction, <I'„(a;) is a random Gaussian variable. 
We see that random variables (j)^^ at different levels 
i describe fluctuations at different length scales. Thus, 
^n{x), containing fluctuations at all length scales, has a 
selfsimilar structure. 




FIG. 1. First two levels of RCS model. 

Next, we calculate the generating function of A^-point 
correlation functions for the RCS model [Eq. (H)]. We 
show that it is quadratic in the sources, i.e., describes a 
Gaussian field theory. We then discuss the relationship 
between the Gaussian field theory for the RCS model and 
the Euchdean field theory C = i(V$)^. 

We consider an arbitrary collection of N points 
large box SI of volume embedded in 
two dimensional Euclidean space. We want to calculate 
the correlation function 



Z„(1,...,7V) =exp 



N 
k=l 



(4) 



where q^, . . . ,qj^ are real and the overline denotes averag- 
ing over all random variables 4>ij ■ The Gaussian random 
variables ^^(a;^,), k = 1, . . . , are not necessarily inde- 
pendent. Their statistical correlations can be measured 

by the scalar product 4; = J27=i[J2'j=i (^^ji^k)'^^ji^l)]■ 
Ij^i measures how "close" Xj. and are in the binary par- 
tition of Q by returning the number /^^ of common boxes 
of volume . . . , /2~^t=i- , respectively, that contain 

both Xj^ and . On the other hand, one can perform 
the Gaussian integrations over all independent Gaussian 



random variables 



and one finds 



Z„(l,...,Af) = e '^S"=i«^ fSi<. 



(5) 



Eq. (|) shows that the RCS model defined by Eq. (|) de- 
scribes a Gaussian field theory, since Eq. (H) implies that 
A^-point functions are constructed from 2-point functions 
in the same way as in the £ = \{V<^f theory. 

Now, let d{xk,xi) = L 2"^'='/^. One can easily show 
that d{xk,Xm) < ina-x[d{xk, xi), d{xi, Xm)], ^Xk,xi,Xm G 
il, i.e., d{x, y) is an ultrametric (and thus a metric), and 
that the 2-point function in the binary tree is given by 



2g , d{x,y) 



GRcs(a;,y) = *„(a;)$„(y) = -|^ln- ^ 



(6) 
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Thus, the RCS 2-point function resembles very much the 
2-point function G{x,y) = In ^^^^^ , in the theory 
^ = 2il(^*)^- Notice that by tuning g = ^f^gA, we 
can match the coefficients in front of the logarithms. 

We have shown that the RCS model yields a Gaus- 
sian field theory, which shares (symbolically) the same 
correlation functions as the C = |^(V<1>)^ field theory. 
The difference between the two theories is that the ultra- 
metric d{x,y) is not the Euclidean metric \x — y\. The 
scaling exponents for the two theories are nevertheless 
the same, despite their different metrics. We believe that 
the multifractal properties of the critical wave function 
for Dirac fcrmions with random vector potential can be 
obtained from the RCS realization. We give support to 
this claim by comparing exact results for the RCS model 
with numerical simulations for the C = 2^i'^^)'^ theory. 

We now turn to the calculation of r(q) within the RCS 
model. It is possible to define box probabilities on our 
random binary tree which obey multifractal scaling of 
the form defined in Eq. (|^). To this end, we define 
local random events at level n in the binary tree by 
On{x) = exp [2<i>„(a;)] [compare with Eq. (||)]. From Eq. 
(0) , we see that 0„ (x) describes a random multiplicative 
process. Be aware, however, that 0„(a;) and Oj^{y), for 
two points X and y, are not independent random vari- 
ables. From On{x) we can construct another random 
variable 



(7) 



where the summation extends over all microscopic boxes 
^Ij of volume = 2~"i^, labelled by any x'^' G flj. In 
analogy to Eq. (|l|) , we then introduce 



lnZ„(g;<l>J ^ lnZ„(l;$J 



(8) 



for any given random event 0„(x) = exp[2<i>„(a;)]. Here, 
we had to divide On{x) by Z„(l; <i>„) to extract properly 
normalized box probabilities: hence the second term in 
the bracket of Eq. (||). Also, the prefactor 2 comes from 
space being two dimensional, whereas the prefactor In 2 
results from having chosen a binary partitioning. As a 
function of q alone, r„(g;$„) shares by construction all 
the analytical properties of multifractal scaling exponents 
for box probabilities We are going to show that the 
large n limit of r„((7;$„) exists and, most importantly, 
is independent of the random event 0,j. In other words, 
T{q) = hm„^oo T„(g; $„) is self averaging. 

To this end, note that we can interpret Z„(g;(f>„) as 
the (random) partition function of directed polymers on 
a binary tree. The inverse temperature is \q\, there are 
2" possible distinct configurations of directed polymers 
which are denoted by all paths (s^, . . . ,s„) on the bi- 
nary tree, and each configuration is weighted by a ran- 
dom Boltzmann weight with random energy 2$„(x[-'l). It 



is clear from Eq. (||) that we are interested in the ran- 
dom variable w„(g;<i>„) = In Z„(q; $„), which is es- 
sentially the (random) free energy per unit length of the 
model of directed polymers in a random medium. The 
asymptotic properties in the limit n — > oo of the random 
variable w„(g; <&„) have been investigated in [ ^Jl5| . 
It is shown in ITq] that almost surely 



lim w„(g;$„ 



In Zi(g;$i) , 



Qc < \l\ 



(9) 



Averaging over the random variables 0^ on the first gen- 
eration of the tree is denoted by an overline. Eq. (|^) 
holds for any disorder such that the annealed average 
Z^{q; $j^) is a well defined even function of q. We have 
restricted ourselves to Gaussian random variables. The 
critical value q^ is defined by the unique minimum of 
q^^ ln[Zj(g; $1)], q>0- The existence of a critical "tem- 
perature" q~^ was anticipated in ||l^. Above the critical 
temperature, the quenched and annealed averages over 
the free energy agree. Below the critical temperature, 
they do not. In Ref. iQ the low temperature phase is 
identified as a glassy phase. 

The Gaussian average on the right hand side of Eq. 
(0) is easily performed. For example 



(10) 



There are two distinct regimes depending on the strength 
g of the disorder. In the weak disorder regime defined by 
q^ > 1, Eqs. (|l^,^,||) yield almost surely 




2(1 -^)2q. 



(2 



9c > kl 



(11) 



in the limit n —>■ 00. In the weak disorder regime, the 
parabolic approximation T*(g) = D*{q){q — 1), where 
D*{q) = 2 — (2/g^)g, is exact for all moments q satisfy- 
ing |<z| < q^- Here, the parabolic approximation (PA) is 
obtained from the annealed disorder average over Gaus- 
sian random variables r* {q) = - [2/ (2 In 2)] [In Z„ (g; <&„) - 
q In Z {1\ <&„)] instead of the quenched average implied by 
Eq. (g) . The PA breaks down for large moments in view 
of the inequality lim„_oo "-~^ln ^„(q; $„) < lnZj(q;<I>J, 
9 1 > 9c — 1- It should be noted that the field theory ap- 
proach of |l^ to Dirac fermions coupling to the disorder 
of Eq. (prelates primary fields with negative scaling di- 
mensions to the multifractal scaling exponents r* {q) (for 
integer valued q) with q^ = 2iT/gj^. In the strong disorder 
regime defined by < 1, the quenched and annealed av- 
erages of Z^[q] <&„) are unequal for all integer moments. 
The PA completely breaks down for integer moments q: 



-2(1 



M), 



q \2 



Ic < 
9c > l9l- 



(12) 
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The multifractal analysis of Eqs. (|11|,|12|) consists in 
performing a Legendre transformation from T{q) to f{a). 
The Legendre transformation is well defined since T(g) is 
concave in the weak as well as strong disorder regimes. 
The interpretation of /(a) is that it yields almost surely 
the Hausdorff dimensions of interwoven Cantor sets char- 
acterizing a typical distribution |l{x,^^yp). Indeed, /(a) 
is strictly positive on its domain of definition, in con- 
trast to the Legendre transform /* (a) of the scaling ex- 
ponents T*(g) (which we recall are calculated from an- 
nealed disorder averages as in [p^). For weak disorder, 
/(a) is a parabola defined on D_^_ < a < D_, where 
Dj_ — 2(1 =F For strong disorder, /(a) is a 

parabola defined on < a < 8/q^. In both cases, the 
parabola takes the maximum value 2, the dimensionality 
of space, and vanishes when a = with finite slopes. 
The finiteness of the slope of /(a) at the end points of its 
domain of definition comes about as a result of the tran- 
sition to linear behavior of r((j) when q> q^. This transi- 
tion is interpreted as a phase transition to a glassy phase 
in the context of directed polymers in random medium. 
This property of /(a) in the RCS construction should 
be contrasted to the deterministic two scale Cantor set 
where the two slopes of f(a) are infinite at Dj^ 

We have compared numerical calculations of the scal- 
ing exponents «((?,$) obtained from the C = ^(V$)^ 
field theory with the exact results derived from the 
RCS construction. We generated four Monte Carlo re- 
alizations for the disorder exp[— J d^a;(V$)^(a;)] on a 
square lattice made of 512 x 512 sites. Within statisti- 
cal fluctuations due to the finite size of the lattice, the 
four Monte Carlo calculations W]vjc(<;, $) and v{q) from 
the RCS model [Eq. (|ll])] agree. 

Finally, we can obtain directly from 2g7(^*^)^ 
the scaling behavior of the sequence of ratios Rm{q)^ 
Z"(q;$)/Z((j;$)"\ Here, Z((j;$)= /c;2a:exp[2g$(x)], 
m is an integer, and averaging is defined by Eq. (H). 
With the help of dimensional analysis, we find that Rm[q) 
contains a term scaling as {a/L)'^"^^'^\ with A,„(<z) = 
2(171 — 1)[1 — ^-^^^f-q'^]. Taking the thermodynamic limit 
yields a finite non- universal number for Rm{q) if 9 < qcrm 
whereas Rm{q) diverges if g > qcm- Thus, there exists 
a decreasing sequence of critical moments = . It 

can be shown along the lines of jl^ that by analytically 
continuing to the limit m = 1, we find the critical 

which is the counterpart of 



.^Dirac 



moment 
Eq. (Il0|). For q 



Dirac gg^j^ commute disorder 



eraging with taking the logarithm of Z{q] and thus 
rely on conventional techniques of field theory to calcu- 
late Tpjj.jj^(g). For q > q^"^'^ , disorder averaging does 
not commute anymore with taking the logarithm, and 
we have to use the RCS model to calculate Tj^-^^^ (q) . No- 
tice that with the identification g^^ = [(47r)/ln2]g, the 
critical moments calculated in the two models coincide. 

All these arguments strongly suggest that Eq. ^ de- 
scribes the exact multifractal properties of the critical 



wave function for Dirac fermions with random vector po- 
tential. Moreover, our results support the conjecture on 
the universality of T(g), (i.e., independence of impurity 
realizations) at criticality in the problem of localization 
in two dimensions j^]. Finally, the methods developed 
in this letter may be just a first manifestation of more 
general principles. We have found a faithful map be- 
tween a particular RCS and a massless free field theory. 
These two theories share symbolically the same correla- 
tion functions, even though they are defined in differ- 
ent metric spaces. Quantities such as scaling exponents 
should not depend on the particular realization of the 
field theory, and thus one has the choice to work with 
whichever is simpler. Whether such maps between field 
theories in different metric spaces can be generalized to 
non-trivial field theories is an interesting open question. 
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